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Asymptotic behavior of singular values of the acoustic 

observation problem 

M. N. Demchenko* 


Abstract 

We consider the problem of recovering of initial data in the IBVP for the wave- 
type equation in the half-space by the solution restricted to the boundary. The 
singular value decomposition of this problem is concerned: the asymptotics of sin¬ 
gular values is obtained. 

Keywords: observation problem, singular value decomposition, spectral asymp¬ 
totics. 


1 Introduction 

Fix an integer d ^ 2 and put Rf = {x = ,Xd ) G R cZ | Xd > 0}. Suppose 

q £ C , Q°(Mi) is a real-valued function. Consider the following initial boundary value 
problem for scalar function u(x,t ), ( x,t ) £ Rf. x R: 

d’fu - Au + qu = 0, d Xd u\ Xd= o = 0, 

u\ t =o = 0, d t u\ t =o = v, (1) 

where v £ CygRf). Introduce the operator 

O : v H- u\? 0 , 

where So = x R (the time-space boundary) and u is the solution of dH) for the 
given initial data v. The operator O is well defined as the solution u is regular. 

Next we introduce a restriction of the operator O. Let U C Rf be a bounded open 
set that satisfies D C Rf. Also let S be a bounded relatively open subset of So, such 
that S C <9R+ x (0, oo). Consider the operator Oq : —> C°°( S) acting as follows 

OEv = (Ov) | E . 

The operator Oq can be continued as compact operator acting from L 2 (D) to .^(S) 
(see sec. ED. 
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We consider the problem of recovering v by OqV (the coefficient q is given), which 
we call the observation problem after JTj. Invertibility of Oq depends on the geometry of 
fl and X. An example of the case when Oq is invertible (more precisely Ker Oq = {0}) 
give E = T x (0,T) and fl satisfying fi C n (U 7e r-S 7 (T)). Here T > 0, T is an 
open bounded subset of cM+, B 1 (T) is an open ball of radius T centered in 7 [ 33 - 
Note that similar problems were also considered in IP, where the method based on 
Carleman estimates was developed. In this paper we are interested in the procedure of 
recovering v based on the singular value decomposition (SVD) for Of r Recall that for 
a compact linear operator A acting in Hilbert spaces Hq —> H\ there exist orthonormal 
basis {v n } C Hq and {/ n } C H i such that 

Av = ^2 s n{A) (v, v n ) Ho f n - 
1 


Here s n (A) are singular values of A defined as s n (A) = (A n (A* A)) 1 / 2 , where A n {A*A) 
are eigenvalues of the compact operator A*A being numbered in non-increasing order 
with the multiplicity taken into account. If Ker A = {0} and / belongs to the range of 
A then 

A~ x f = ^s n (A) _1 (f,f n ) Hl v n . (2) 

71^1 

It can be seen that the stability of inversion (in particular) depends on the behavior 
of s n (A). In this paper we prove that s u (Oq) decrease as a negative power of n (re¬ 
lation (gD). Note that we do not need the invertibility of Oq and thus no geometric 
conditions on Q and X are imposed. 

To study the asymptotics of s n (OE) we approximate the operator Oq by Fourier 
integral operators (FIO), which yields approximation of (Oq)*Oq by pseudo differential 
operators (TDO) of order —2. Then we apply the result on asymptotics of singular 
values for non-elliptic TDOs of negative order [5]. The theory of FIOs is a natural tool 
for analysis of operator Oq. Furthermore, FIOs were applied to the problems of integral 
geometry (8j, 0, 15( CEZl [20] (the list of references is not complete), which are similar to 
the problem of inverting Oq. 

The principal symbols of 'kDOs that approximate (0^)*0^ vanish on the subset of 
T*Q \ 0 of nonzero measure. If ( Oq)*Oq was precisely an elliptic TDO (i.e. operator 
with nonvanishing principal symbol) then ((Oq)*Oq)~ 1 would be a TDO of order 2; 
if additionally Ker Oq = {0} then the operator (O^) -1 would be bounded in certain 
Sobolev spaces since 

(oEr 1 = ((oEyoEr^oEy. 

In our case the operators that approximate (0^)*0^ do not have even parametrices, 
so the observation problem is ill-posed. The principal symbols of approximating TDOs 
vanish outside of the conic set AZ C T*H\0 which is defined in sec.[2j The same situation 
arises in problems of integral geometry with limited data - see [aEunuHEo]. The 
conic set AZ is called an audible zone (or visible zone ) there. The results of these papers 
show that it is possible to reconstruct some part of singularities of unknown function, 
namely - the intersection of the wavefront set of the function with the audible zone. 
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This gives a hope that series ([2]) applied to Oq allows to reconstruct v microlocally in 
the audible zone, the reconstruction being stable in some reasonable sense. However, in 
the present paper we only obtain the asymptotics of s u (Oq). 

Note that for some problems of integral geometry singular values and singular basis 
(i.e. v n and f n in ([2]) ) were found in a nearly explicit form. It concerns some particular 
(though very important) cases such as limited angle tomography [ 12 L [13] and exterior 
Radon transform [IS]. It seems that there are no such results for Oq even for q = 0 and 
some certain P and X. 

Acknowledgements. The research was supported by RFBR grants 14-01-31388 
rnol-a, 15-31-20600 mol-a-ved, 14-01-00535 and SPbGU 6.38.670.2013. 
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2 The audible zone and the asymptotics of s n (0 q) 

To define the audible zone AZ and formulate the main result we introduce some notation. 
A point in the plane So will be identified by (a/, t), x' € M“ -1 , t £ R (we use identification 
of (M+ and M rf_1 ). For x £ R d we put x' = (xi,...,a^-i). We also use natural 
identification of cotangent bundle T*P and Introduce the following mapping 

from T*P to Xo: 

i(y,v) = (y'-yd^,Vdj^ , y e P, y e T y *P 

(the function 7 is defined if rjd 7 ^ 0, i.e. almost everywhere in T* O). Note that the 
first component of 7 (y, 77 ) is an intersection point of the line {y + sy, s € M} and dW+. 
while the second component is the distance between this intersection point and y. The 
audible zone is defined as follows 

AZ = {(y, 7 ) € T*P \ 0 | Vd + 0, 7 (y, 7 ) G X}. 

We say that K C T*f2 \ 0 is a conic set if it is invariant with respect to the mapping 
(y, rj) 1 —x (y, rr)) for any r > 0. The set AZ is conic, and besides, it is invariant with 
respect to the mapping (y,rj) i-x (y, — rj ). 

Denote by n(A, A) the number of singular values s n (^4) greater than A > 0. The 
main result of the paper is the following Theorem. 

Theorem 1. Let P C M+, X C cM°|_ x (0, 00 ) be open bounded subsets satisfying P C M+, 
X C cM+ x (0, 00 ), /id(<9X) = 0 (pd is the d-dimensional Lebesgue measure on Xo, 9X 
is taken in the topology 0 /X 0 ). Then 

limA'yog.A) = (2 *r d JJ A /yw (^ -') • P) 

6 is the Heaviside function. 
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Note that the integral in d3J is finite. Indeed, the Heaviside function vanishes if 7] 
is large and \rjd\ ^ 1; if r/ is large and \r]d\ < 1 then (y,rj) ^ AZ since y G H and H is 
separated from 8M.+ . We conclude that the Heaviside function is nonzero only on the 
bounded set of ( y,r /). 

The result ([3j implies the following behavior of singular values: 

Sn(O^) ~ (0 ' , n -X oo, (4) 

where cr is the r.h.s. of ([3]). Note that the asymptotics dU has sense if a ^ 0. Also 
relation (J3J) may be written in the equivalent form in terms of (Oq)*Oq: 

]im\ d / 2 n«OEyoE,\) = a. (5) 

A—>0 

Note that here n coincides with the counting function of (positive) eigenvalues of (Oq)*Oq, 
i.e. the number of eigenvalues greater than A. The relation (|5j) is a Weyl-type asymp¬ 
totics for positive operator (Oq)*Oq. 

In the rest of the paper we prove Theorem |T] 


3 Operators I, I± 

Here we introduce and investigate FIOs connected with the problem dTJ) . Some basic 
concepts and facts of the theory of FIOs are used (canonical relations, composition of 
FIOs, principal symbol of composition of FIO with its adjoint) - we send the reader 
to books mmm for details. First we express Ov in terms of the following Cauchy 
problem for u(x,t) in l^xl (i.e. in the whole space): 

d 2 u — A u + qu = 0, 

u\t=o = 0, d t u\ t =o = v. (6) 

Here q(x) = q(x f , |xd|) is an even continuation of q to the whole space (q is smooth as 
q is compactly supported in M^), and v € Cq^M^). Suppose v £ Co°(R^) and take v 
such that v = v in and v = 0 in R fi \ R^_. Then the solutions of (JT]) and (J6|) are 
related as follows 


u(x, t) = u(x', Xd, t) + u(x' , — Xd, t), Xd > 0, t € R. (7) 

Indeed, u(x' , —Xd,t) satisfies the wave-type equation in dTJ for Xd > 0 and has zero 
Cauchy data for t = 0 in R^_. The boundary condition for Xd = 0 is obviously satisfied. 
It follows from dZJ that 

Ov = 2u\x 0 . ( 8 ) 

Next we represent the solution of ([SJ in terms of FIOs: 

u = I + v + I-V (modulo a smoothing operator), 
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(I±v)(x,t) = (2tt) d [[ dydy e t<f>± ^ x,t,y ’ v ^ a±(x, t, y, y) v(y), 

JJ R d xE d 

where <p±(x,t,y,y ) = (x — y)rj ± t|r/| is a nondegenerate phase function with phase 
variable y and a± is an amplitude of order —1. The amplitude a± is classical: 


d± ~ 

J >1 


; 0 ‘) 


( 9 ) 


£ is a smooth function such that £(s) = 1 for s > 2 and £(s) = 0 for s < 1 , functions 
a±) are homogeneous of degree —j in y. Functions are solutions of certain transport 
equations m, but we need only a± : 


4 1} = H^\v\) 


-i 


( 10 ) 


Note that the definition of a FIO in [ 2 T j implies the normalization factor (27r) _rf_1 / 4 
rather than (27r)~ d in formula for I±. Nevertheless our choice will be more convenient 
in further considerations. 

The operators I± are continuous from to C' 0 O (R d+1 ) (as well as from £'(R d ) 

to 2T(R rf+1 )). Hence we may consider the restriction I±v |s 0 for v € Co°(R d ). The 
restriction can be written in the form of FIO as well: 

(I±v)(x',t) = ( 2 vr) _d ff dydrje lip± ( x '’ t ’ y ’ r,) a±(x',t,y,ri)v(y), (x',t)<£S 0 , 

JJ K d xK d 

where ip± and a± are just the restrictions of (p± and a±. Note that 

V±(x',t,y,y) = x'rf - yy ± t\y\ 


is a nondegenerate phase function. The mapping 

v (->• (I+v + i-v - it)|s 0 = I+v + i-v - ii|v;o 
is smoothing. In view of ((HJ) we obtain the following representation of O : 

O = 2 (1+ + /_)) (modulo a smoothing operator). (11) 

Now let x € Co°(£o) be a real-valued function, and put 

I± = | (yoo Qgd ) • (12) 

We consider I± as a FIO mapping functions in R^ to functions in Ho- The phase function 
ip± and the amplitude x a ± of the operator I± are functions of (x',t,y,y), where (x',t) 
and y vary over So an< l R d \ {0} correspondingly (just as for I±), while y varies over 
R^_. 

For the operator / = 2 (1+ + /_) due to (fTTT) and (fT2l) we have 

\0 = I (modulo a smoothing operator). (13) 
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In sec. [5] some approximation of the characteristic function of the set £ will be chosen 
as the multiplier % in (1121) . This will provide that the restriction of the operator xO 
to functions supported in hi will approximate the operator Oq, which means that the 
analogous restriction of I will approximate Oq. 

The canonical relation of I± looks as follows 

C± = {{x\ t, dj't'pi, y , -d y <p±) | dr,<p± = 0} C (T*£o \ 0) x (: T*R d + \ 0). 

Note that —d y (p± = r/, hence 

C± = {(x',t,d x > tt ip±,y,r]) | d v (p± = 0 } 

and the phase variable rj may be treated as an element of T*W^_. Introduce an open 
conic subset of T*M+ \ 0: 

K = {{y,v) G T*R d + \0\r, d ± 0}. 

Consider the pair of canonical transformations from K to T*£o \ 0: 

(y,v) •-t (x?,t,£), (xf,t)='y±(y,7i),g = r/,£d = ±\'nl ( 14 ) 

where 

7±(y>*7) = (V -yd.—, ±yd —) • 

V Vd m J 

Analysis of ip± shows that the canonical relation C± is a graph of the mapping (I14|) . 

Due to the multiplier x i n the definition (1121) we may consider the compositions I±I± 
and I^I±. Since C± is a canonical graph the composition I±I± is a TDO (see [21]): 

(I±I±v)(y) = (2tt )~ d [[ dydy e ^ y ~ y ^ b±(y,rj)v(y). 

J 

Since I± and I± are classical FIOs (i.e. Schwartz kernels can be represented by oscillatory 
integrals with classical amplitudes) the composition I±I± is a classical \kDO: 

&±~5^C(M)&±\ (15) 

where b± are homogeneous of degree —j in y. The adjoint operator I± is a FIO with a 
canonical relation C±. Composition Itl± is a smoothing operator in since CLoC± = 
0 . 

Now we may write 

I*I = 4(/+/+ + III-) (modulo a smoothing operator) 
and I*I is a classical TDO in of order —2 with symbol 


6 = 4 ( 6 + + 6 _). 


( 16 ) 
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From da it follows that 


(xO)*xO = I*I (modulo a smoothing operator). (17) 

(t£\ 

We now calculate the principal symbol of the 'I'DO I±I Z t, he. the function b±(y,r]) 
in (1151) . We use the formula from [211 sec. 8.6] to express the principal symbol in terms 
of ip± and a±. Consider the map 

H- {y, ~d y (p±{x',t,y,ri),d ri (p±(x',t,y,rff) = (y, rj, d v tp±(x', t, y, rj)) 


and denote by A± the absolute value of its Jacobian determinant. Then the principal 
symbol of I±I± in (y. rj) E K equals 


(I°± ) I 2 a± 1 )(7 ±(y,v),y,v), 


where is the leading term of the amplitude a±, due to (flOj) we have d± = ±(2i\rj\) 1 
Direct calculation shows that for (y,rj) E K we have A± = |?7d|/|?7|, so 


b±\y, rj) 


x(7±(y,v)) 2 

Mv\ \vd\ 


The principal symbol b'± should be continued with zero to (T*M^ \0) \ K. Due to (fl6l) 
for the principal symbol b l 2 ' in the expansion 


j> 2 


(18) 


we have 

x(n +(y, v)) 2 + x(i-(y, v)) 2 

\v\\m\ ' 1 j 

Further we will apply operators O, I to functions in L, 2 (£l), which is possible since 
these operators can be expressed in terms of FIOs that act on compactly supported 
distributions in Ml and D is separated from 3M+. 



4 Spectral asymptotics of ^DOs 

We use the following result on spectral asymptotics of a TDO with symbol C(\y\)d m ^ (y. rj) 

(Bv)(y) = (2vr) _d f f dydy e l{y ~ y ^ C(|y|) 6 (m) (y, rj)v{y), 

J Jflx R d 

where is homogeneous of degree —m (m > 0) in and b E x (M d \ {0})), 

function ( is the same as in Q. B is a compact operator in L 2 (D). The result of [5] 
claims that 


lim A d/m n(5, A) 

A—>0 


(2 tt) 1 [[ dydr)d(\d m) (y,ri)\ - 1). 

J JflxR d 


( 20 ) 
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Note that B is not supposed to be self-adjoint or elliptic in [5], Also note that [5j 
concerns the more general case when is anisotropically homogeneous and is only 
continuous in 12 x (R d \ {0}). 

Now let /q be the restriction of I to L2(12). The result (1201) implies the following 
asymptotics for I^In (the operator I^In is compact since it is the restriction of the 
TDO I*I of order -2 to L 2 (S2)) 

limA d / 2 n(/a/n,A) = (27r)- d ff dydy 9(\b^(y,y)\ - 1), (21) 

JJnxR d 

where b^ is given by (1191) . To prove (1211) we need to show that 

lim A d / 2 n(I^In, A) = lim A d ^ 2 n(B, A), (22) 

A—>o A->0 

where B is a TDO in 12 with symbol C{\ r l\)b^{y, ff)- The relation (1221) claims merely that 
terms b^\ j ^ 3, in the expansion (fl8l) do not influence the asymptotics of n(I^lQ, X). 
Although this fact seems to be trivial, the author could not find an appropriate reference, 
so a short proof is provided here. We establish the estimate 

A d/3 n(I^I n -B,X)^C VA > 0, (23) 

which means that n(I^Ia — B, A) is estimated by less power of 1/A than (1/A) d / 2 and (1221) 
will then follow (see [4] sec. 11.6]). 

Choose an open bounded set 12 7 containing 12, such that 12' C R+, <912 7 £ C°°, and 
denote by Iqi the restriction of I to L2(12 7 ). Let B' be the TDO in 12' with symbol 
({\v\)b {2) (y, rf). Denote by Pq the projector onto L 2 (Q) acting in L2(12 7 ). We have 

n(I* n I n - B, A) = n(P n (i£,/ n , - B’)P n , A). (24) 

Choose p £ Co°(f2'), p\n = 1, and put F = p(I^,Iq/ — B')p. Since 

P<i(I*ykv - B')P a = P n FPn 


we have 

n(Pn(It :v kv - B')P n , A) < n(F, A). (25) 

Here we used the following inequality 

Sn(AL), s n (LA) ||L|| s n (A) 

for a compact operator A and a bounded operator L. Put D = (1— A) 3 / 2 , where A is 
the self-adjoint Laplace operator in Q 7 with Dirichlet boundary condition. D is a TDO 
in 12 7 [19]. The composition FD is also well-defined as a TDO of order ^ 0, hence it 
is a bounded operator in ^(D 7 ). The operator D~ l is also bounded in L2(12 7 ). For 
v £ L 2 ( 12 7 ) we have 

(Fv,Fv)l 2{Q ') = {FDD- l v,FDD~ l v) L2(QI) ^ WFD^^D^v^)^). 
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Hence F*F < \\FD\\l 2{nl) D- 2 and 

n(F*F, A) ^ n(\\FD\\l 2ini) D- 2 , A) = n(D~ 2 , X/\\FD\\ 2 L2{ni) ). 

The well known asymptotics of eigenvalues of —A yields 

n(D~~ 2 , A) ^ CX~ d ^. 

We arrive at 

n(F, A) ^ CX~ d/3 . 

In view of f[Mj) . (1251) this leads to (1231) . 

5 Proof of Theorem |1] 

In this section we prove the relation ([5]). 

First show that OE is a compact operator from L 2 (D) to L^iTX). Since the operator 
IqIq is compact (see sec.|4]) the operator Iq acting from ^(D) to Z/ 2 (So) is also compact. 
Denote by Oq the restriction of O to L 2 (D). Now due to (fT3l) the composition x@£i is 
a compact operator, which implies that Oq is also compact since we can choose x such 
that x|e = 1 . 

To prove ([5]) we will make different choices of function x i n the definition (|12l) . First 
suppose that 0 ^ x ^ 1, x = 1 on E and suppx C cM+ x (0, oo). For v € L, 2 (£l) we 
have 

(Oqv,0%v ) l2 ( S ) ^ (xOnv, x Onv) L2{ So) . 

Therefore (C>E)*OE ^ (xOnYxOn and so 

n((OE)*OE, A) < n((xOnTxOn, A). (26) 

It follows that 

liin A rf / 2 rj((0fj)*O^, A) < Um 'X d / 2 n(( X O n YxOn, A). (27) 

A— >o A— >o 

Due to (fl7l) the “principal parts” of operators (xOnYx^ii and coincide, hence 

lim A d/2 n((xOn)*xOn, A) = lim A d/2 n(I^In, A) = cr 2 , (28) 

A—>0 A—>0 

where <72 is the r.h.s. of (1211) . To prove this one should repeat the proof of (1221) . 

Now put 

n{y,v) = x{i+{y^)) 2 + x(i-(y,'n)) 2 , (y,v) = xxh(y,y)), 

Xs is a characteristic function of £. We have 

< k < 1. (29) 
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Indeed, the first inequality is obvious if 7(y, 7) ^ E. Suppose 7(7, 7) £ E, then 7 + (y, 7) £ 
E in case rjd > 0 and 7_ (y, 7 ) £ E in case r]d < 0. Now the inequality ks ^ k follows from 
x|s = 1- The second inequality follows from the fact that conditions 7 + (7,7) £ suppy 
and 7_(y,r/) £ suppy can not hold simultaneously since suppy C <9K°|, x (0, 00). 

Due to the first inequality in (1291) and formula (fT9l) we have a ^ 07 . Now we need 
to estimate 02 — o. 

Since D C there exists e > 0 such that if y £ Q and |%|/| 7 | < e then n(y, rj) = 0. 
If |r/| > e -1 / 2 and y £ D then n(y, v)/(\ 1 l\\ r ]d\) < 1- Indeed, if |/?d|/| 7 | ^ £ (otherwise 
ft(y, rj) = 0) then due to (1291) 


«(y,7) 1 

I 7 ?! \Vd\ ^ I 7 ?! M < 

This means that the integrals in (1211) and in ([3]) are taken over the set y £ D, \rj\ < e -1 / 2 . 
The difference <72 — <7 can be estimated by the measure of the set of (y, 7) such that 

M < e” 1/2 > Ks(y,ri) < rz{y, rj. (30) 

Due to (1291) the second condition implies that ky, ( 7 , 7 ) = 0, and so 7 ( 7 , 7 ) ^ E, which 
means that neither 7 + (y, 7 ) nor 7 _(y, 7 ) belongs to E. However, on the set (130j) we have 
/c(y, 7 ) > 0, and so 7 +(y, 7 ) or 7 _(y, 7 ) belongs to suppy. We conclude that every (y, 7 ) 
from the set (1301) satisfies 

| 7 | < £ -1/2 , 7 + ( 7 /, 7 ?) or 7 -(2/) r /) belongs to suppy \ E. (31) 

Now choosing y such that supp y shrinks to E (note that in this case £ does not depend 
on y) we make the set (1311) shrink to the set of ( 7 , 7 ) such that I 7 I < £~ 1//2 and 7 + (y, 7 ) 
or 7 _(y, 7 ) belongs to 3E. This set has zero measure in H x W l since fid(dT,) = 0. 
This means that choosing appropriate y we can provide that the difference 07 — a is 
arbitrarily small. Together with (1271) and (|28p this means that 

Ihn A d / 2 n((0^)*0Q, A) ^ a. (32) 

A—>0 

Taking y £ C'q°(E), 0 ^ y ^ 1, we obtain the inequality reverse to PUI) . therefore 
we have the relation (which is a counterpart of (1271) 1 

lim A d / 2 n((0f))*0Q, A) ^ lim A d / 2 n((yOn)*yO n , A). 

A-s-0 A-s-0 

Then arguing the same way as in proof of (1321) we obtain that 

lim A d/2 n((0%)*0%, A) ^ a. 

A—>0 

Together with (l32l) this yields ((5|) and thus Theorem |T] is proved. 

Now we expose the reason why we used the smooth multiplier y in the definition (1121) 
of the operator I Instead of this we could try to deal with I' n := 2(/ + + /_) considered 
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as an operator from ^(O) to L 2 (E) (see formula (Hill ). It seems to be possible to 
consider (J^)*!^ as a 'I'DO with the principal symbol 

As (7 (g. 7 )) 

M M 

(xs is the characteristic function of £), which is discontinuous in 0 x (R (i \{0}). However, 
in this case the result (1201) of [5] can not be applied since the principal symbol is required 
to be continuous there. Note that classical results on spectral asymptotics of \HDOs with 
discontinuous symbols by H. Widorn and their improvements also can not be applied to 
Mn) *Iq - we do not go into details here. 
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